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ABSTRACT 

We analyse the cluster color-magnitude relation (CMR) for early-type galaxies in two of the richer 
clusters in the z ~ 0.9 supercluster system to derive average ages and formation redshifts for the 
early-type galaxy population. Both clusters were observed with the Advanced Camera for Surveys 
aboard the Hubble Space Telescope through the F606W and F814W filters, which brackets the rest- 
frame 4000 A break at the cluster redshifts of z ~ 0.9. We fit the zeropoint and slope of the red cluster 
sequence, and model the scatter about this relation to estimate average galaxy ages and formation 
redshifts. We find intrinsic scatters of 0.038— 0.053 mag in {Veoe—Isu) for the E and E-I-SO populations, 
corresponding to average ages of 3.5 — 3.7 Gyr and formation redshifts Zf = 2.4 — 2.6. We find at least 
one significant difference between the C11604-I-4304 and 4321 early-type CMRs. 4321, the less X-ray 
luminous and massive of the two, lacks bright L* ellipticals. We combine the galaxy samples to fit a 
composite CMR down to 0.15L*, and find that the slope of the combined cluster CMR is significantly 
steeper than for RX J0152. 7-1357 but consistent with MS 1054-03, both at similar redshift. The slope 
of the C11604 CMR at the bright end (L > 0.5 L*) is fiatter and consistent with the CMR slopes found 
for other high redshift clusters. We find evidence for increasing scatter with increasing magnitude 
along the early-type CMR, consistent with a 'downsizing' scenario, indicating younger mean ages with 
decreasing galaxy mass. 

Subject headings: galaxies: clusters: individual (CL1604) - galaxies: elliptical and lenticular - galaxies: 
evolution 



1. INTRODUCTION 

The majority of elliptical galaxies appear to have 
formed the bulk of their stars at very early times. 
In the local universe, elliptical galaxies over a 
wide range of environment follow a tight color- 
mag nitude rela tion (CMR) with ext remely smal l scat- 
ter dy isvanathan & Sandaec 1977; Bo wer et al.l 119921: 
iHogget al.-2004:.McIntosh et al.-2n051 In a galaxy clus- 
ter, there is an obvious bimodality in galaxy colors that is 
well-correlated with galaxy morphology. The 'blue cloud' 
is dominated by spiral and irregular galaxies, and the 
prominent ridge of red galaxies is composed mainly of 
ellipticals and SOs. More massive red sequence galaxies 
are on average older, more metal-rich, and contain rela- 
tively more products of core-collapse supernovae, indicat- 
ing more rapid formation, than lower mass red s equence 
galaxies (jBernardi et al .^ 2005: Nelan et a l. 2005). 

The early-type cluster CMR, or red cluster sequence, 
is such a regular feature that it can be used to find rich 
galaxy clusters at z > I (e.g. Gladders & Yee 2005). The 
origin of the slope of the CMR is generally accepted as 
a mass-metallicity relation, where more massive galax- 
ies have higher metallicities (e.g. Kodama & Arimoto 
1997). The scatter about the CMR is related to age 

^ Department of Physics and Astronomy, Johns Hopkins Uni- 
versity, 3400 North Charles Street, Baltimore, MD 21218. 

^ Department of Physics and Astronomy, Washington State Uni- 
versity, Pullman, WA, 99164-2814. 

3 STScI, 3700 San Martin Drive, Baltimore, MD 21218. 

** UCO/Lick Observatory, University of California, Santa Cruz, 
CA 95064. 

"" European Southern Observatory, Karl-Schwarzschild-Strasse 2, 
D-85748 Garching, Germany. 

^ Leiden Observatory, Postbus 9513, 2300 RA Leiden, Nether- 
lands. 



differences at a given galaxy mass. The tightness of the 
CMR and its lack of observed evolution with redshift im- 
plies that CMR galaxies formed at high redshift, z > 2. 
Since galaxies change color more quickly at younger ages, 
pushing CMR observations to higher redshifts allows us 
to more accurately determine the formation redshifts of 
early- type galaxies. 

The ACS Intermediate Redshift Cluster Survey (Ford 
et al. 2004) covers 8 clusters in the redshift range 0.8 < 
z < 1.3 with the Advanced Camera for Surveys (ACS; 
Ford et al. 1998) aboard the Hubble Space Telescope as 
part of the ACS Guaranteed Time Observations (GTO) 
program. Of the previous results, the most relevant to 
the current study are the analyses of the CMR in the z = 
1.24 cluster RDCS J1252. 9-2927 (Blakeslee et al.l2003bD 
the z = 1.1 cluster RDCS J091 0-^5422 IjMei et aL I 
l2006aj) . t wo clusters in the z = 1.27 Lynx supercluster 
CL0848 l)Mei et al. Il2006bj) . and two z ~ 0.83 clusters , 
MS 1054-03 and RX J0152. 7-1357 fBlakeslee et al."2006). 
These four studies spanning 6 clusters found small intrin- 
sic scatters of 0.03 — 0.06 mag about the CMR, average 
luminosity- weighted ages of 2.2 — 3.5 Gyr, and forma- 
tion redshifts of z > 2.3 — 2.7 in a concordance ACDM 
cosmology. 

In this paper we extend these early-type CMR 
studies to include two clusters at the low-mass 
end of previously studied clusters. We present an 
analysis of two low X-ray luminosity clusters in 
the CL1604 supercluster, C11604-H4304 and 4321, 
which have been previously studie d in so me detail 
iCunn et all 119861: iCastander et all Il994t fOke et alJ 



1998t iPostman et alJ 119981: iStanford et alJ Il998t' 



Lubin et a,lJl2nnnt IPostma,n et al.1 120011 IStauford et alJ 
20021 IHolderi et a,lJ I2004L Wide field imaging of 



the surrounding regions uncovered red galaxies with 
colors consistent w ith passive galaxies at z ~ 0.9 
IjLubin et al.1 120001) . and follow-up spect roscopy con- 
firme d the existence of a supercluster l|Gal fc LubinI 
Hooi)- C11604-h4304 has a velocity dispersion of 
962 ± 141 km s^^ from 67 redshifts l|Gal fc LubinI 
|2004|) . and a relatively low temperature and lumi- 
nosity: Tx = 2.51+j^-^^ keV and Lx,boi = 2.0 x lO^" 
ergs s-1 l|Lubin et all |2004|) . C11604- H4321 has a 
veloc ity dispersion of 640 ± 71 km s~^ (|Gal fc LubinI 
linni- C11604-h4321 was not detected by ROSAT, 
and has an upper limit on the 0.1 — 2.4 keV X-ray 
luminosity of L^ < 4.76 x lO'*'^ erg s~^, compared 
to the 0.1 - 2.4 keV detection of C 11604-H4304 of 
L^ = 7.43 ± 1.59 X lO''^ erg s'^ IjPostman et al.l 
120011) . Based on their measured velocity dispersions, 



C11604-f 4321 is at least a factor of ~ 2 less massive than 
C11604-I-4304. 

The paper is organized as follows. In § 2 we describe 
the observations, object selection, and photometry. In § 3 
we present our color-magnitude diagrams. § 4 discusses 
the results of the CMR fits, § 5 and 6 discuss the CMR 
scatter as a function of magnitude and comparisons to 
other work, and in § 7 we present our conclusions. 

2. OBSERVATIONS AND REDUCTIONS 

CL1604-f4304 and CL1604-h4321 were observed with 
the ACS Wide Field Channel as part of a guaranteed 
time observation program (proposal 9290). These clus- 
ters were observed for 2 orbits in each Vgoe and I^u 
at a single pointing. Th e data were pro cessed with 
the Apsis pipeline (Blake slee et al.ll2003aj) . Our pho- 
tometry is calibrated to the AB magnitude system us- 
ing zeropoints from Sirianni et al. (2005). Object de- 
tection an d photometry was performed by SExtractor 
(jBertin fc _ Arnouts 1996 ) incorporated within the Apsis 
pipeline. A more d etailed description can be found in 
iBenitez et aP l|200^ . We use MAG.AUTO when quot- 
ing total broad-band magnitudes, a s this is the best esti - 
mate of a galaxy's total magnitude (" Benitez et alJl20M) . 
Zeropoints for the Vqqq and /814 filters are 26.486 and 
25.937 AB magnitudes, respectively. We correct the 
galaxy magnitudes and colors for Galactic extinction: 
this amounts to 0.033 ± 0.005 in Fgoe and 0.021 ± 0.003 
in Isii for CL1604+4304, 0.039 ± 0.006 in Veoe and 
0.024 ± 0.004 in Igu for CL1604+4321. 

Extensive vis ual morpholog y catalogs were cre- 
ated by MP l|Postman et all I2005D . Morpholo- 
[ :ies were determined visua lly on the T-type system 
de Vaucouleurs et al.lll99E . All galaxies in the field 
with /814 < 24 magnitude were classified. Approxi- 
mately 10% of the galaxies were also classified by three 
independent classifiers to estimate the classification er- 
rors. Majority agreement was achieved for 75% of ob- 
jects with Z775 < 23.5. There was no significant sys- 
tematic offset between the mean classification from the 
in dependent class i fiers. More information can be found 
in lPostman et all l|200,5f) . A T-type -5 < T < -3 corre- 
sponds to elliptical galaxies, — 2 < T < — 1 corresponds 
to lenticular (SO) galaxies, T=0 to SO/a galaxies, and 
T > to Sa and later- type galaxies. 

We adopt Af^ = -21.3 AB mag fror n the Schecter 
funct ion fit to galaxies in the 2dF survey IjNorberg et al.l 
12002(1 and -1.04 magnitudes of luminosity evolution at 



z = 0-90 ((van Dokkum fc Stanford! l2003t Hfolden et alJ 
12005ft . This corresponds to an observed magnitude of 
/g;^4 = 21.9 at z = 0.90 for an early-type spectrum 
with Vgoe ~ -^814 = 1-8, using the conversion /gi4 = 
Mb + 43.820 + 0.089 x {Veoe - hu) - 0.794 obtained 
from the BC03 solar metallicity single burst models with 
a Salpeter IMF. Thus, our visual morphological classifi- 
cation extends to 0.15 L* . Throughout this paper we use 
Ho = 70 km s-i Mpc~\ n^ = 0.30, and J7a = 0.70. 

We select all galaxies with T < brighter than /gi4 = 
24.0 mag. Because early-type galaxies have color gra- 
dients that can bias the measurement of the CMR, we 
measure galaxy colors within fixed spatial apertures de- 
fined in the reddest filter. We produce 80 x 80 pixel 
cutouts and subtract the background image produced by 
Sextractor. Other objects in the cutout are masked, and 
we fit a single psf- convolved Sersic profile with GALFIT 
(Peng et al.ll2002() . with the position, orientation, ellip- 
ticity, and effective radius as free parameters. The n 
index was constrained to be between 0.1 < n < 4, and 
the 'sky' was fixed at 0. We 'clean' the cutouts with the 
CLEAN algorithm to remove the blur ring of the PSF 
l(Hogboml I197I iBlakeslee et al.l l2003bl) . The residuals 
are added back to the 'cleaned' images to conserve flux, 
and photometry is performed within a circular aperture 
of rh = Re\/b/a, derived from fitting the F814W image. 
If the derived rh is less than 3 pixels, we set it to 3 pixels. 
Color errors are derived from the RMS images generated 
in the Apsis pipeline (Blakeslee et al. 2003a). These im- 
ages have estimated errors per pixel, including Poisson 
flux errors, read noise, cosmic ray rejection, dark, and 
bias frame errors. The RMS errors should be close to 
but slightly smaller than the true errors, due to corre- 
lated errors as a result of the geometric correction. 

3. THE COLOR MAGNITUDE RELATION 
We fit linear relations of the form: 



(Veoe - hii) = slopeilsu ~ 23. 



Co 



where cq is the zeropoint and slope is the slope. For 
each cluster we fit the full early-type sample (E, SO, and 
Sa), as well as the Es, the SOs, and the Es-|-SOs. For 
C11604-I-4304 and the combined cluster sample, the fits 
were restricted to galaxies with colors 1.5 < T^oe ^-^814 < 
2.0, then used an iterative linear fit with 2.5cr clip- 
ping based on Bisquare weighting (Press et al. 199^) 
and est imated the scatte r with a biweight scale esti- 
mator ((Beers et al.lll990() . We note that linear least- 
squares fitting with 2.5cr clipping produced similar re- 
sults. The 4321 CMR was less straightforward to fit than 
the C11604+4304 CMR. Restricting the inital sample to 
1-5 < ^^606-414 < 2.0, 1.6 < V"606-4i4 < 2.0, or within 
0.06 X 3.0 mag of the C11604-h4304 CMR aU lead to dif- 
ferent CMR slopes and scatters. The incompleteness of 
our spectroscopic sample for this field hampers our abil- 
ity to accurately measure the 4321 CMR. In the end we 
chose to perform an initial clip around the C11604+4304 
CMR, then used an iterative linear fit with 2.4(7 clipping 
based on Bisquare weighting. However we discuss CMR 
results obtained from other initial clipping choices. 

Uncertainties in the slopes, zeropoints, and scatters 
were derived with 10,000 bootstrap simulations. To esti- 
mate the intrinsic scatter about the fitted CMR we sub- 
tracted, in quadrature, the average galaxy color error 



from the measured scatter. We also determined the mean 
scatter required to make the reduced ^ of the hnear fit 
equal to 1. This method yielded intrinsic scatters lower 
by « 10 — 20 %, and are the values we quote. 

The scatter about the CMR can be used to constrain 
the formation epoch of the CMR galaxies. Using a simple 
model that assumes that the difference in age is relatively 
constant as a function of magnitude (galaxy mass), the 
intrinsic scatter can be mod eled in three ways. Follow - 
ing the method first used by Ivan Dokkum et all l)200Cl() , 
the first method assumes galaxies form in single bursts 
at random times between t^iz — 5) and tend- As we in- 
crease iend, the mean age of the 10000 model galaxies 
increases, and the distribution of galaxy ages also nar- 
rows. We assume a minimum ie„d of 0.5 Gyr. Galax- 
ies need at least 0.5 Gyr after ceasing star formation to 
redden and arrive on the red cluster sequence. These 
models assume no correlation in galaxy formation times 
within a cluster halo. The relationship between intrin- 
sic scatter, luminosity- weighted mean galaxy age, and 
iend is shown in Fig. |21 for models with a Salpeter IMF 
and three metalhcities (Z=0.008, 0.02, 0.05; Bruzual & 
Chariot (2003)). Another method assumes galaxies form 
stars at constant rates from random times t\ and t^ be- 
tween io and iend- Again we vary tend from 0.5 Gyr 
from the time of observation to to- A third method is 
to assume exponentially declining star formation rates, 
starting and ending at random times t\ and ti between 
to and tend. Previous studies have found that the ages 
derived from this constant star formation model are sys- 
tematically lower by 10 - 20% (Blakeslce ct al. 2003b). 
This is approximately what we find here for intrinsic scat- 
ters of 0.05 — 0.03 mag. However the largest changes in 
mean ages would be from metalhcities lower than solar, 
according to the Bruzual & Chariot 2003 models. As 
is known from high-resolution spectroscopy in the local 
universe, massive elliptical galaxies have solar or super- 
solar metalhcities, although it is not clear how smoothly 
the metallicity varies with luminosity along the red se- 
quence. These models would predict a decreasing scatter 
with luminosity if the mean ages did not depend on lu- 
minosity. 

In the following sections we quote constraints on galaxy 
ages and formation times from only the single burst 
method with solar metallicity. Errors on the ages and 
formation redshifts come from adding and subtracting 
the bootstrapped Icr errors to the intrinsic scatter to de- 
rive the asymmetric lu limits on the average ages and 
formation redshifts. 

We compare our fitted CMR relations to the relation 
for Coma from Bower et al. (1992). We transformed the 
l] — V slope and zeropoint to U ~ B using the BC03 so- 
lar and super-solar 12-13 Gyr SSP models. From this 
we derived the relation: S(u - B) = 0.66A(f/ - V). 
Using the empiric al Kinney -Calzetti elliptical, SO, and 
Sa templates (Kin nev et al.1. 1996) and the Coleman el- 
liptical template ( Col eman et al.l 1198(1) gave a similar 
result. We transformed the U — B slope, —0.054, to 
the Veoe — I&u slope of —0.066 at z ~ 0.90 with the 
derived relation ^(T/eoe - hu) = l-23(5(t/ - B). This 
transformation depends somewhat on the templates used 
to derive the colors. With the empirical templates we 
derived the relations 5{U - B) = 0.51A(f/ - V) and 



We convert the observed Vgoe — -^814 colors measured 
in AB magnitudes to rest-frame U-B measured in Vega 
magnitudes using the IRAF SYNPHOT package and 
the BC2003 Z=0.008, 0.02, 0.05 metallicity models with 
single bursts and ages between 1 — 6 Gyr. We per- 
formed a linear fit to these models and find 6{U — B) = 
0.79A(T/606 - hu) at z = 0.9. 

3.1. Redshift Color Corrections 

The color difference in Vgoe ~ Isu between a redshift 
of z = 0.90 and a redshift oi z — 0.92 depends strongly 
on the underlying galaxy SED, if one considers the full 
range of galaxy types. For example, a galaxy matching a 
z = SO template spectrum would be 0.0041 mag redder 
sX z — 0.92, while an Sa template spectrum would be 
0.0130 mag redder, but an Sb template would be 0.006 
mag bluer. The Veoe— -^814 colors of the CMR galaxies fall 
between the Sa (Vgoe — ^8i4 = 1.82 at z = 0.90) and Sb 
(^606 —^814 = 1-51 at z = 0.90) empirical templates, but 
the sh ape of the SED is no t likely to match. Instead we 
use the lBruzual fc Charloti l|2003) single bursts and short 
exponentially declining SFR models and we find correc- 
tions of the order of A(y606 - hu) = 0.02 - 0.8Az. The 
range in redshift for each cluster is Az ~ 0.01, or max- 
imum corrections of 0.008 in (Vgoe — -^814) color. These 
corrections are small enough that we neglect them. 

3.2. CLI6O4+43O4 

The color-magnitude diagram of C11604+4304 is 
shown in Figure |3^. The red cluster sequence is promi- 
nent in the color-magnitude diagram of this cluster field. 
The fitting results are summarized in Table ^ For the 
E-t-SO sample we find 

(VeoG-hu) = (-0.076±0.009)x(/8i4-23)-l-1.738±0.008 

and an intrinsic scatter of 0.038 ±0.004 mag. The slopes, 
zeropoints, and scatters of the E-|-SO-|-Sa, E-f-SO, E, and 
SO relations are consistent within the errors. 

In Figure 01 we show the color magnitude diagram 
with elliptical galaxies as circles, SOs as squares, SO/a as 
stars, and late-type galaxies as triangles. Spectroscop- 
ically confirmed cluster members are indicated by open 
boxes. Two of the brightest four ellipticals are spectro- 
scopically confirmed cluster members. Galaxies included 
in the CMR fit are indicated as larger filled symbols. 
The CMR relation is overplotted as a solid line, and the 
transformed Coma relation as a dashed line. Modeling 
the scatter about the red cluster sequence for the E-I-SO 
sample with the stochastic single bursts, we find an av- 
erage galaxy age of 3.711q'o3 Gyr, with an average for- 
mation redshift z/ = 2.60J:°;°i^. The Es and SOs do not 
have significantly different intrinsic CMR scatters, so we 
find similar ages and formation redshifts. 

3.3. CLI6O4+432I 

Visually, there are many more blue galaxies in the 
CL1604-F4321 field than in the C11604-h4304 field (see 
Fig-EJj the red cluster sequence appears less defined, and 
it lacks bright, /814 < 21.5 magnitude elliptical galaxies. 
The CMR fitting results arc summarized in Table El Us- 
ing an inital clip of 0.06 x 3 mag around the C11604-f 4304 
CMR, we find for the E-hSO sample 
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^814J 



(-0.045±0.021) X (/8i4-23) + 1.782±0.015 



and an intrinsic scatter of 0.053 ± 0.010 mag. For the 
E-only sample we find a slope of —0.066 ± 0.027 and an 
intrinsic scatter of 0.048 ± 0.011 mag. Again for this 
cluster, the slopes, zeropoints, and scatters of the E-only 
and SO-only relations are the same within la errors. 

For the E+SO intrinsic scatter, with the single burst 
models we find an average age of S.SIq'j Gyr, and Zf = 
For the elliptical-only sample of 19 galaxies. 
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we find an average age of 3.60_Q]^g Gyr, and a mean 

Zf — Z..O_Q^Q. 

If we use the same color cuts and sigma-clipping 
when fitting the C11604-f4304 and 4321 CMRs, the 
4321 fitted CMR slopes are significantly flatter than the 
C11604-h4304 CMR, and the intrinsic scatter is signif- 
icantly larger than the C11604+4304 scatter. For ex- 
ample, using 1.5 < Vgoe — hu < 2.0 and 2.4(T clip- 
ping, the E-I-SO sample the slope and intrinsic scatter 
are 0.00 ± 0.04 and 0.066 ± 0.009. Shown in Fig. El this 
CMR includes elliptical galaxies redward of the Coma 
slope, and SOs with Veoe ~ -^814 — 1-5 — 1.6. Some or all 
of these galaxies may be interlopers. On the other hand, 
they may be cluster members which have yet to evolve 
to a tight CMR, and we are biasing the CMR slope and 
underestimating the intrinsic scatter with our initial cut 
about the C11604-I-4304 CMR. The E-only CMR slope 
and intrinsic scatter are -0.022±0.037 and 0.062±0.011, 
still significantly flatter and with more scatter than the 
E-only slope initially clipped around the C11604-I-4304 
CMR. 

If we try to cut based on distance to reduce contamina- 
tion by non-cluster members, we still find flatter slopes. 
If we define the center of the cluster as the location of the 
brightest elliptical (16:04:33.61, -^43:21:04.0, J2000) and 
fit the E-I-SO CMR within 0.5r2oo, we find a slope and 
intrinsic scatter of -0.034 ±0.022 and 0.043 ±0.010. But 
this spatial cut excludes a spectroscopically confirmed el- 
liptical galaxy (z=0.9274). Going out to 0.7r2oo includes 
this galaxy, and we flnd a slope and intrinsic scatter of 
-0.019 ±0.028 and 0.055 ±0.015 (we exclude 1 SO galaxy 
with Veoe — ^814 = 1.59). We list these different possi- 
ble fits for the 4321 CMR to highfight the uncertainty of 
CMR fitting when the number of CMR galaxies is small, 
and spectroscopic membership is unknown. In this case 
the initial color and magnitude cuts have a large impact 
on the derived CMR. 

3.4. Combined Cluster CMR 

We can combine the galaxy samples because of the 
small redshift difference between the clusters. The dif- 
ference in V(306 — ^814 color between an early-type galaxy 
at z = 0.90 and z — 0.92 is small enough to be negligible 
(0.0008 AB mag for BC03 tau models). The fitting re- 
sults are summarized in Table |21 The combined cluster 
CMR is shown in FigureEl with the E+SO+Sa fit as the 
solid line, and the Coma relation as the dashed line. For 
the combined cluster E+SO sample we find 
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606— -'814J 



(-0.068±0.010)x(/8i4-23.) + 1.754±0.009 



and an intrinsic scatter of 0.052 ± 0.006 mag. The com- 
bined CMR and this fit is shown in Figure [7| For the 
E-only sample, we find 

(^606-/814) = (-0.070±0.013)x(/8i4-23.) + 1.759±0.011 



and an intrinsic scatter of 0.052 ±0.007 mag. We find the 
same slopes, zeropoints, and scatters within the errors for 
all galaxy samples. 

From the intrinsic scatter about the CMR for the E+SO 
sample including 65 galaxies, we find an average age of 
3.50lo;|!)7 Gyr, and 2/ — 2.4lQ;Qg with the stochastic sin- 
gle burst model. The E-only sample of 42 galaxies has 
an intrinsic scatter of 0.052 ± 0.007 mag, identical to the 
E+SO sample, so we also derive t = 3.5 Gyr and Zf — 2.4. 

The combined cluster CMR has many spectroscop- 
ically confirmed cluster members, including 8 of the 
brightest 14 CMR galaxies. Fainter than /814 ~ 22.2 
we note that the spectroscopic members tend to be bluer 
than the CMR, likely due to the combination of color and 
i?— band selection. There are only three SO/a galaxies in 
the fit, and they are all fainter than /814 = 23.5. The 
CMR appears to be populated fairly evenly by E and 
SO galaxies in magnitude; no significant segregation is 
seen. As for the individual clusters, the measured CMR 
slope, zeropoint, and scatter for the E-only and SO-only 
combined populations are consistent within the errors. 

4. CLUSTER COMPARISON 

In this section we compare the slopes, zeropoints, and 
scatters of the C11604+4304 and 4321 CMR relations. 
As stated in the previous section, we can directly com- 
pare the CMRs because of the small redshift difference 
between the clusters. 

For the E+SO sample, we find slopes consistent at the 
- 1.4cr level (-0.076 ± 0.009 vs. -0.045 ± 0.021), with 
the formally steeper slope belonging to C11604+4304. As 
discussed in § 3.3, the 4321 CMR is less defined and more 
problematic to fit. A difference in CMR slopes is not 
ruled out, given the incompleteness of the spectroscopic 
samples. For the elliptical galaxies only, the slopes are in 
agreement: -0.074+0.012 vs. -0.066±0.027), as are the 
SO-only slopes (-0.080 ± 0.016 vs. -0.045 ± 0.042). 

The slope of the CMR can be used as a diagnostic 
of the average age of the galaxies IjKodama fc Arimotol 
119971 IGladders et al.llT99l . Within 4 Gyr of formation, 
more metal-rich galaxies become redder faster than more 
metal-poor ones, and an initally flatter slope becomes 
steeper. Thus, if the CMR galaxies formed indepen- 
dently of mass (magnitude), within 4 Gyr of formation, 
the slope should be flatter than its 2; = value, because 
the brighter galaxies are bluer with respect to their final 
colors than the less luminous galaxies. 

To get an idea of the range of slopes one would expect 
at this redshift, in Figure El we show the C11604+4304 
CMR with our fitted CMR relation as a solid line, the 
4321 relation as a dash-dot line, and three models from 
Kodama & Arimoto l|Kodama fc Arimotcll997) with for- 
mation redshifts of Zf = 2, 3, and 5. The model CMR 
relations are summarized in Table 01 The slopes be- 
come steeper as the formation redshift increases, mostly 
from the bright end of the CMR becoming redder. All 
three model slopes are too flat to match the faint end of 
the C11604+4304 CMR relation, but the z/ = 2 model 
matches the CMR relation of 4321. If the zeropoint of 
4321 were not redder than C11604+4304, then we would 
interpret the marginally flatter slope as a consequence 
of a larger formation redshift for 4321. Basically, the 
slope for 4321 is flatter because the galaxies are redder 
at a given magnitude, which is inconsistent with younger 



ages. 

The zeropoints of the E+SO and E-only relations are 
significantly different between C11604+4304 and 432f 
(2.1-2.6cr). The Clf 604+4304 zeropoint is bluer in both 
cases, by about 0.04 magnitudes. This can be seen in Fig- 
ure This is caused by faint galaxies in C11604+4321 
that are redder than the faint galaxies in C11604+4304 
((Ko6 --fsw) = 1-7-1.9 vs. 1.5-1.7). If they are clus- 
ter members, this would imply they are either older or 
more metal-rich than galaxies of corresponding bright- 
ness in C11604-H4304. Given C11604-h4321's lower mass 
and lack of bright ellipticals, it would be surprising if the 
low-mass end of the CMR contains older or more metal- 
rich elliptical galaxies. We think it more likely that the 
low mass end of the C11604+4321 is not yet or only par- 
tially formed, and the faint, red, elliptical galaxies are 
interlopers, pushing the zeropoint redder, and the slope 
flatter. More extensive spectroscopy of the C11604 su- 
percluster region indicates that there is a background 
structure at z = 1.1 - 1.2 in the C11604+4321 field (R. 
Gal, private communication). 

The intrinsic scatters for all samples are not signifi- 
cantly different between the two clusters, implying sim- 
ilar ages and formation redshifts. This is the result we 
derive if we initially clip the 4321 sample around the 
C11604-I-4304 CMR. If we instead use a similar color cut 
and 2.4(7 clipping, then 4321, the lower mass cluster, has 
a much larger scatter, corresponding to an average age 
difference of - 0.4 Gyr between C11604-I-4304 and 4321, 
for both the E-I-SO and elliptical-only populations. 

5. CMR SCATTER AS A FUNCTION OF MAGNITUDE 

With all the evidence that less massive galaxies form 
their stars at later time s than more mass ive galaxies, 
the 'downsizing' concept IjCowie et al.lll99(j) . assuming a 
constant scatter about the CMR as a function of magni- 
tude may not be correct. If the downsizing phenomenon 
also occurs for CMR galaxies, then not only should the 
average color become bluer with magnitude (decreasing 
mass), but the range in colors should also increase, in- 
creasing the measured scatter with magnitude. We inves- 
tigate the dependence of the CMR on galaxy luminosity 
(mass) by first fitting the CMR in restricted magnitude 
ranges, then fixing the CMR slope and zeropoint from the 
fit to the bright end and calculating the intrinsic scatter 
about this relation. 

There are too few galaxies in the C11604-I-4321 
sample for meaningful measurements, so we focus 
on C11604-t-4304 and the combined sample. For 
CL1604-^4304, the intrinsic scatter about the CMR for 
galaxies with 20.5 < /814 < 22.5 (bright end; Mb ~ 
-22.5 to -20.8) is 0.025 ±0.005 magnitudes, and 0.039 ± 
0.005 mag for the range 22.5 < Isi4 < 24. (faint end; 
Mb 20.8 to -19.5), a 2a difference. The slope dif- 
ferences are suggestive, but not statistically significant: 
-0.052 ± 0.045 at the bright end, and -0.096 ± 0.020 at 
the faint end. 

For the combined cluster sample, at the bright end, we 
find 

iVeoe-hu) = (-0.040±0.024)x(/8i4-23.)-hl.780±0.029 
and an intrinsic scatter of 0.030 ± 0.005 mag. There are 
no so/a galaxies in the bright sample. 

We also fit the faint end and restricted the initial sam- 
ple to be between 1.5 < (Veoe — -^814 < 1.95). For the 



E+SO sample, we find a steeper slope of —0.089 ± 0.024, 
a zeropoint of 1.759 ± 0.010, and an intrinsic scatter of 
0.045 ± 0.006. The elliptical-only scatter also increases 
with magnitude, from 0.024 + 0.004 mag to 0.043 + 0.009. 
These differences in scatter are significant at the ^ 2a 
level. 

For the combined cluster population, we fix the slope 
at —0.068 and zeropoint at 1.754, as we found for the 
combined E+SO population, and measure the scatter at 
the bright and faint end about this relation. For the 
E+SO population, we find intrinsic scatters of 0.034 + 
0.005 mag at the bright end, and 0.048 + 0.006 mag at the 
faint end, a 1.8cr result. For the elliptical-only population 
fixing the slope and zeropoint at -0.070 and 1.759, we find 
0.030 + 0.005 mag at the bright end, and 0.050 + 0.007 
mag at the faint end, a 2.3cr result. 

Summarizing these results, with both methods, 1) an 
independent fitting of the slope, zeropoint, and scatter of 
the CMR at bright and faint ends, and 2) fixing the slope 
and the zeropoint from the full fit and compute the scat- 
ter about this relation for the bright and the faint ends, 
we detect marginally significant {2a) increases in scatter 
with magnitude in the E+SO and elliptical populations. 
The difference in scatter amounts to an age difference of 
~ 0.15 Gyr with the solar metallicity single burst models. 

Similar re sults have been previously reported. 
iTanaka et alJ (J20O4) found increasing scatter about the 
CMR found in groups at z — 0.83, although they 
did not find evidence for this in their observed cluster, 
RXJ0152. 7-1357, which we speculate is due to large pho- 
tometric errors relative to ACS. A dependence of CMR 
scatter on CMR galaxy p roperties has been f o und i n the 
studies of MS 1054-03 bv lvan Dokkum et all l|200nl) and 
iBlakeslee et alJ 1)20061) . Restricting the sample to early- 
type galaxies brighter than ^775 < 22.5 magnitude and 
excluding merger candidates, the rest-frame U — B scat- 
ter is 0.024+0.006 mag. But including merger candidates 
and all early-types down to 1775 = 23 magnitudes, they 
found an intrinsic scatter of 0.052 + 0.009 mag. 

An increase of scatter and steepening of the slope com- 
pared to 2: = would be expected if fainter galaxies are 
on average younger. If all CMR galaxies formed at the 
same time, but with metallicity varying with mass, there 
should be no increase in scatter, or possibly a decrease 
in scatter. An increase in scatter at faint magnitudes 
is expected if lower mass galaxies form their stars at 
later times, and the expected offset from the CMR of 
the higher mass galaxies should be bluewards. 

6. DISCUSSION 

C11604+4304 and 4321 are at the low-mass end of pre- 
viously studied clusters at z ^ 1, so the comparison of 
the early-type CMR to other high redshift clusters should 
give us important information about the variation in for- 
mation times of cluster early-type galaxies. A significant 
uncertainty in comparing the CMR slopes between clus- 
ters is in the transformation from the observed slope to 
rest-frame \5{U — B)/5B\, as this is dependent on the 
spectral templates used for the transformation. For con- 
sistency with our comparison sample, we use the same 
templates used for the C11604 slope and scatter transfor- 
mation (see §3) to transform the observed CMR slopes 
for th e elliptical-only sampl es for MS 1054 a nd CI 0152- 
1357 lIBlakesleeet al.ll200l . C10910+5442 lIMei et al. I 
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I2006a|) C11252-292 7 (iBlake slee et"alll20031j) . and Lynx 
E+W l|Meiet al. 1 120061^1. We list the slopes and the 
transformations in Table [S] We note that the quoted 
slope measurements are made over different absolute 
magnitude ranges, which may lead to different derived 
CMR slopes. Another complication is in the filter trans- 
formations, i.e. how well is the underlying galaxy spec- 
tral energy distribution matched by the templates used 
to derive these transformations. 

The three clusters at z - 0.9, MS 1054-03, RX J0152.7- 
1357, and C11604, span a large range of slopes, and 
it is unclear if this is significant. The MS 1054 
and RX J0152. 7-1357 slopes were fitted with a magni- 
tude limit about 1 magnitude fainter than C11604, and 
were fitted using only spectroscopically confirmed clus- 
ter members. If we fit the C11604 elliptical-only sam- 
ple at the same limiting magnitude, we find a much 
shallower CMR slope {\S{U - B)/5B\ = 0.027 ± 0.014) 
that is in agreement with both RX J0152. 7-1357 and 
MS 1054-03. For C11604, the fitted CMR slope at the 
bright end is c onsis tent wit h other clus ters at in terme- 
diate redshift l|van Dokkum et_aU|2000|: [Blakeslc e et al.l 
120031^ lM;r^rarn200fia.kTR]a,keslee et a,1..200fil^. but 
the galaxies at the faint end are relatively bluer, push- 
ing the composite slope steeper. The only cluster with 
a comparably steep slope is Lynx W, however, only 6 
galaxies contribute to this slope. Trends of CMR slope 
with cluster mass and X-ray properties will be discussed 
in Mei et al. (2006c). 

The average formation redshifts of early-type galax- 
ies derived here, Zf = 2.4 — 2.6, agree well with those 
from previous cluster studies: Zf ~ 2.2 a t z = 0.83 
llBlakeslee et al.ll2006j) . Zf « 2.3 at z = 1.10 llMei et al I 
I2006a|^ . and~~ 2.8 at z = 1.27 l|Mei et al. ll2006bD . 
We measured the scatter about the CMR in restricted 
magnitude ranges, and find that scatter increases as a 
function of ma gnitud e. 

iWake et al.Tl)2005|) studied 12 clusters at z - 0.3 cov- 
ering a factor of 100 in X-ray luminosity, and put an 
upper limit of 2 Gyr on the average variation in CMR 
galaxy ages from the lack of observed dependence of 
CMR zeropoint on X-ray luminosity. Our results are not 
inconsistent with this, as an age range of 2 Gyr is quite 
large, much larger than the range of < 1 Gyr that we 
probe here. They also noted that the 5 least X-ray lumi- 
nous clusters (0.1 — 0.7 x lO^^erg s~^) had sparse CMRs, 
forcing them to combine clusters to measure a composite 
faint X-ray cluster CMR. We find a similar qualitative 
difference in the color- magnitude diagrams and CMRs 
of C11604-H4304 and 4321. There is a lack of bright, 
Mb < —20.9 elliptical galaxies in 4321, possibly because 
mass assembly/BCG formation in 4321 is incomplete. 

The CMR scatters and slopes of the z — 1.27 Lynx 
E and W clusters were found to depend on distance 
from the cluster center l)Mei et al. Il2006bri . increasing 
and steepening as the distance increased. We find no 
such dependence for either cluster sample, but the num- 
ber of galaxies is small and the radial coverage is slightly 
less than 1 Mpc. We also find no difference in the 
CMR relations for the E and SO galaxies separately. 
So far, only one study has found a differ ence in the 
E and SO cluster populations at z ~ 1. IjMei et al. I 
l2006ai) found that the CMR for SOs in C10910 (z = 1.10, 



Lx.Boi = 2.5i!^j X lO*"* erg s^^; Stanford et al. 2002b) 
had a significantly bluer zeropoint than for the ellipticals, 
implying younger ages or lower metallicities at similar lu- 
minosity. 

The intrinsic scatter for the combined cluster sample 
agrees well with other clusters at similar redshifts, con- 
firming a lack of measured evolution of scatter about the 
CMR down to low cluster X-ray luminosities. This lack 
of evolution in the measured sca tter could be due to a 
comb ination of 'progenitor bias' ijvan Dokkum fc Franxl 
|2001|) and the fitting technique. We found that the in- 
trinsic scatter increases as a function of magnitude, and is 
significantly larger for CMR galaxies from 0.15 — 0.4L* in 
our combined cluster sample. The combination of galax- 
ies within 0.5— 2 L* having similar properties over a range 
of environment, and fewer CMR galaxies with L< 0.4L* 
present at z ~ 1 (e.g. De Lucia et al. 2004) could 
'conspire' to kee p the measured sca tters constant. As 
also discussed in iHolden et alJ (|2£)04) , by sigma-clipping 
around the CMR, we are selecting the oldest cluster 
galaxies at any given redshift, not the cluster population 
that will evolve into the z = CMR population. 

Although the CMR is o bserved in large area galaxy 
surve ys out to z = 1 (Be U et alJ 120041: iMcIntosh et alJ 
I2005J) . we know that the CMR evolves. Striking evidence 
for the formation of the CMR and its depe ndence on en- 
vironment was recently demonstrated by iTanaka et alJ 
l)2004j) . While there is a clear CMR in the field, groups, 
and clusters at z = 0, the field CMR is not present 
at z = 0.83, and becomes less defined in groups at 
z = 0.83 and m > M* + 2. This was detected as an 
increase in CMR scatter as a function of increasing mag- 
nitude. Building on this theme, we find no evidence 
for evolution, other than passive, in the bright end of 
the CMR (L > 0. 5L*), in agreement with many othej 
studies (van Dokkum et alJl200 0: Blakeslc e et al. 2003^^ 
Holden et al.ii2004l: JMei et al. i r2006a b: Blakeslee et al 



|2Q0(|). However, considering CMR galaxies down to 
0.15L*, we find tentative evidence for a dependence of 
CMR slope on cluster X-ray luminosity or mass, in the 
sense that less luminous clusters have steeper slopes. If 
this is confirmed by further studies, it must be due to a 
combination of 'downsizing', where the fainter, less mas- 
sive CMR galaxies form later than the brighter, more 
massive ones in a given cluster, and an overall younger 
age for faint CMR galaxies in the lower mass clusters. If 
'downsizing' were not invoked, meaning all CMR galaxies 
form simultaneously, then younger clusters should have 
fiatter slopes. 

Earlier studies of the scatter in lower redshift clus- 
ters could only put a lower limit of Zf > 1 — 1.5 for 
the average formation redshift of the CMR galaxies (e.g. 
van Dokkum et al. 1998). Now that we can observe 
galaxy clusters in this redshift range, we can put much 
stronger constraints on their formation redshifts. The 
range z = 3 — 5 is emerging as crucial for understanding 
the formation of massive cluster ellipticals and their tran- 
sition from rapidly star-forming to passive. Although 
semianalytic and hydrodynamic models of galaxy for- 
mation are now able to re produce a c olor bimoda lity of 
galaxies from z = — 1 IjDave et al . 2005: Kang et all 
120051 iMenci et all 120051 they still produce bright blue 
galaxies at z = 0, the SFRs of most model red galaxies 
decline only slowly after z = 2, and is driven by major 



merging events l|Menci et al.ll2005|) . The star formation 
decline is much too slow to reproduce the tight CMRs 
found in z ~ 1 clusters, and major merging including gas 
(dissipative mergers) is not observed to occur for a sig- 
nificant fraction of massive cluster galaxies. What these 
models probably do get right is that the color bimodality 
of galaxies is driven by major differences in their progen- 
itor star formation histories at z = 4 — 5, which we can 
hopefully test by continuing to push our observations 
outward in redshift. 

7. CONCLUSIONS 

We presented an analysis of the color-magnitude rela- 
tion in C11604-H4304 (z = 0.90) and C11604+4321 (z 



0.92V both part of the 16h supercluster l|Gal fc LubinI 
12004ft . Based on their X-ray luminosities, both of these 
clusters are at the low mass end of the range of previously 
studied clusters at z = 0.8 — 1.3. Our main conclusions 
are as follows. 

• Both C11604-h4304 and 4321 have identifiable 
early-type color-magnitude relations. However, 
C11604-I-4321, the less massive of the two clusters, 
lacks bright M* ellipticals, and we suggest that it 
has not yet had time to assemble bright early-type 
red sequence members. From our analysis of the 
CMR relations, we find a redder zeropoint for 4321. 
The CMR slope and scatter depend strongly on our 
choice of initial color cuts, and the possibility of 
significant differences in slope and scatter between 
C11604-f4304 and 4321 is open. 

• The slope of the combined cluster CMR is signifi- 
cantly steeper than for RX J0152. 7-1357 but con- 



sistent with MS 1054-03, both at similar redshift. 
This may be due to the difference in magnitude 
limit used for the fitting, the different observed fil- 
ters, or a real difference in the galaxy populations. 
The slope of the CMR at the bright end (0.5-2 L*) 
is flatter and consistent with the CMR slopes found 
for other high redshift clusters. 

• From the intrinsic scatter about the CMR for the 
full sample of 65 E-f SO galaxies, we find an average 
galaxy age of 3.50^0 J? Gyr, and an average z/ — 
2 4+°-^ 



• We find evidence for an increasing intrinsic scatter 
with increasing magnitude, previously only known 
in the group and field environment at these red- 
shifts. This implies an decreasing mean galaxy age 
with decreasing galaxy mass, of about 0.15 Gyr 
from - IL* to - 0.3L*. 
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TABLE 1 
CL1604+4304 






Sample 


Ngal 


Zeropoint 


Slope (Tint 


5{U - B)^/SB 


<T^nt{U-B), 


E+SO 

E 

SO 


38 
23 
15 


1.738±0.008 
1.739±0.011 
1.738±0.012 


-0.076±0.009 0.038±0.004 
-0.074±0.012 0.038±0.006 
-0.080±0.016 0.033±0.007 


-0.061±0.007 
-0.059±0.009 
-0.064±0.012 


0.030±0.003 
0.030±0.004 
0.026±0.005 



Note. — For galaxies brighter than /814 — 24 mag. Errors on tlie rest-frame U-B slopes and scatters 
do not include errors from the transformation. 



TABLE 2 
CL1604+4321 



Sample ^gal Zeropoint 



Slope 



S(U-B)JSB a,„t{U-B), 



E+SO+Sa 


30 


1.769+0.015 


-0.053+0.021 


0.059+0.008 


-0.042+0.015 


0.043+0.006 


E+SO 


27 


1.782+0.015 


-0.045+0.021 


0.053+0.010 


-0.036+0.015 


0.040+0.007 


E 


19 


1.783+0.018 


-0.066+0.027 


0.048+0.011 


-0.053+0.020 


0.036+0.009 


SO 


9 


1.47+0.63 


-0.045+0.042 


0.052+0.024 


-0.036+0.031 


0.039+0.018 



Note. — For galaxies brighter than /814 — 24 mag 
do not include errors from the transformation. 



Errors on the rest-frame U-B slopes and scatters 



TABLE 3 
Combined Cluster Sample 



Sample ^^gal Zeropoint 



Slope 



5{U-B),/5B a^„t{U-B), 



E+SOH 
E+SO 
E 
SO 



Sa 



69 


1.753+0.009 


-0.065+0.011 


0.055+0.006 


-0.052+0.008 


0.044+0.005 


65 


1.754+0.009 


-0.068+0.010 


0.052+0.006 


-0.054+0.008 


0.042+0.005 


42 


1.759+0.011 


-0.070+0.013 


0.052+0.007 


-0.056+0.010 


0.042+0.006 


22 


1.740+0.010 


-0.070+0.013 


0.041+0.006 


-0.056+0.010 


0.033+0.005 



Note. — For galaxies brighter than /si4 — 24 mag. For the SO fit we used an initial color cut of 
1.55 < Vaoe - /8i4 < 1.95. 



TABLE 4 
Kodama-Arimoto Models 



z formation CMR Slope CMR Zeropoint 



-0.036 
-0.047 
-0.047 



1.78 

1.81 

1.84 



TABLE 5 
Rest-frame CMR elliptical-only slopes 



Cluster 


Redshift 


Observed Slope 




transformation 




—S{U - B)/5B— 


Reference 


Coma 




-0.079 + 0.007 


slope{U - 


-B) = 0.66 slope{U - 


V) 


0.052 + 0.005 




(1) 


RX J0152.7-1357 


0.83 


-0.007 + 0.011 


slope{U - 


- S) = 1.19 slope{re25 


-«775) 


0.008 + 0.013 




(2) 


MS 1054-03 


0.83 


-0.042 + 0.014 


slope{U - 


-B) = 0.91 slopeiVeofi 


-J775) 


0.038 + 0.013 




(2) 


CI 1604 


0.91 


-0.070 + 0.013 


slope{U - 


-B) = 0.79 slope(Vfiofi 


-/814) 


0.055 + 0.010 






CI 0910 


1.10 


-0.033 + 0.015 


slope{U - 


- B) = 1.12 slope{i775 


- Z850) 


0.037 + 0.017 




(3) 


CI 1252-2927 


1.24 


-0.020 + 0.009 


slope{U - 


- B) = 1.72 slope{irrs 


— 2:850) 


0.034 + 0.015 




(4) 


Lynx E 


1.27 


-0.025 + 0.012 


slope{U - 


- B) = 1.79 slope{i775 


- Z850) 


0.045 + 0.021 




(5) 


Lynx W 


1.27 


-0.043 + 0.015 


slope{U - 


' B) = 1.79 slope{irr5 


- 2850) 


0.077 + 0.027 




(5) 


Note. — (l)IBower et al.lIT993). (2)|Blakeslcc ct al. (2006:), 


(3) Mei ot al.jk2006ail. 


f4)iBlakcslee ct al. 12003b), 


(5)11 


ei et al. 1 
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Fig. 1. — Rcst-framc \5{U — B)/5B\ CMR slopes vs. cluster redshift. 
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Fig. 2. — Intrinsic scatter in Veoe — ^814 at z = 0.9 vs. luminosity-weighted mean galaxy age and tend with different models. From left 
to right: (a) single bursts with Z=0.008 (black), Z=0.02 (medium gray), Z=0.05 (light gray), (b) exponentially declining star formation 
rates with r = 1 and solar metallicity, (c) constant star formation rates between ti and t2 chosen at random between to and t^„j. In all 
three models, the minimum tg„^ is 0.5 Gyr (z=1.02). 
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Fig. 3. — Color magnitude diagrams of the cluster fields: C11604+4304 (left) and C11604+4321 (right). 
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Fig. 4. — CL1604+4304 (left) and CL1604+4321 (right) color magnitude relations. Visual morphological types are indicated by circles (Es) 
and squares (SOs). Red sequence members are indicated by red circles (Es) and orange squares (SOs). Open boxes indicate spectroscopically 
confirmed cluster members. Open circles indicate spectroscopically confirmed interlopers. The small dots are early-type galaxies that are 
not included in the fit. These include the galaxies rejected by the iterative fit and the spectroscopically confirmed interlopers. The CMR 
fit to the E+SO sample is overplotted as a solid line. We do not plot the other (E, SO, E+SO+Sa) relations because they are so similar to 
this relation. The black dashed line is the transformed Coma relation. 
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Fig. 5. — CL1604+4321 CMR, using the same color cuts and sigma clipping as for C11604+4304, instead of an initial cut around the 
C11604+4304 CMR. Symbols are the same as for Fig. g] The solid line is the fit to the E+SO sample. 
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Fig. 6. — CL1604+4304 with three CMR models from Kodama & Arimoto with formation redshifts of 2, 3, and 5 (dashed lines, from 
bluer to redder). Galaxy symbols are the same as in Figure The solid lines are the C11604+4304 and 4321 E+SO CMRs (+4304 is 
steeper). The black dashed line is the transformed Coma relation. 
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Fig. 7. — The combined cluster sample. Galaxy symbols are the same as in Figure|^ including SO/a galaxies as orange stars. The CMR 
fit to the full sample is overplotted as the solid line. Ths fit to only the Es is the dotted line. The dashed line is the transformed Coma 
relation. 



